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Abstract. Tliis article reports the first optical frequency measurement of the IS — 3S transition in hydrogen. 
The excitation of this transition occurs at a wavelength of 205 nm which is obtained with two frequency 
doubling stages of a titanium sapphire laser at 820 nm. Its frequency is measured with an optical frequency 
comb. The second-order Doppler effect is evaluated from the observation of the motional Stark effect due 
to a transverse magnetic field perpendicular to the atomic beam. The measured value of the lSi/2(i^ = 
l)-3Si/2(i^ = 1) frequency splitting is 2 922 742 936.729 (13)MHz with a relative uncertainty of 4.5x 10"^^ 
After the measurement of the IS — 2S frequency, this result is the most precise of the optical frequencies 
in hydrogen. 



1 introduction 

The hydrogen atom has a central position in the history of 
atomic physics. As it is the simplest of atoms, it has played 
a key role in testing fundamental theories, and hydrogen 
spectroscopy is associated with successive advances in the 
understanding of the atomic structure. Since the advent in 
the seventies of tunable lasers and methods of Doppler free 
spectroscopy, hydrogen spectroscopy has been renewed in 
the last decades. Consequently, several optical frequencies 
of hydrogen are now known with a fractional accuracy 
better than 10~^^. In a long series of experiments, Hansch 
and coworkers have improved the precision on the mea- 
surement of the IS — 2S frequency to obtain now a rel- 
ative uncertainty of about 1.4 x 10^^^ |T|. In Paris, we 
have studied in the nineties the 2S — nS/D two-photon 
transitions with n = 8 and 12 (n is the principal quantum 
number) [2I3I4J . For instance we have measured the fre- 
quency of the 2Si/2~8D5/2 transition with an uncertainty 
of 5.9 kHz, i.e. a relative uncertainty of 7.6 x 10~^^. The 
goal of these high precision measurements is to determine 
the Rydberg constant i?oo and the hydrogen Lamb shifts. 
The hydrogen energy levels can be conventionally ex- 
pressed as the sum of three terms: the energy given by the 
Dirac equation for a particle with the reduced mass, the 
first relativistic correction due to the recoil of the proton 
and the Lamb shift. The first two terms have an exact 
expression as a function of the quantum numbers and of 
the fundamental constants (the Rydberg constant Hao , the 
fine structure constant a and the electron-to-proton mass 
ratio rrie/mp). The Lamb shift takes into account all the 
other corrections: corrections due to quantum electrody- 
namics (QED), other corrections due to the recoil of the 



proton and the effect of the proton charge distribution. 
The calculation of the Lamb shift is very difhcult and a re- 
view of the results obtained so far is made in the report of 
the CODATA (Committee on Data for Science and Tech- 
nology) [5]. Today, for the IS level, the uncertainties are 
due to the calculation of the two- loop and three- loop QED 
corrections (this uncertainty is estimated to 3.7 kHz) and 
to the measurement of the proton charge distribution. Us- 
ing the value of the radius rp of this charge distribution 
deduced from the electron-proton scattering experiments 
(rp = 0.895 (18) fm) ^, this uncertainty is 50 kHz. At this 
level, the uncertainty of the theoretical value of the Lamb 
shift is not limited by the ones of R^o, o: and rrie/mp. 
Consequently, the uncertainty in the proton radius rp is 
a strong limitation to extract the Rydberg constant from 
high precision measurements. For instance, using the fre- 
quency of the IS — 2S transition and the calculated value 
of the IS and 2S Lamb shifts, the Rydberg constant can 
be deduced with a relative uncertainty of only 1.8 x 10~^^ 

In practice, it is possible to avoid this difficulty by us- 
ing the 1/n^ scaling law for the Lamb shift. Numerous 
terms of the Lamb shift vary with the principal quan- 
tum number exactly as 1/n^ (for instance the effect of the 
charge distribution of the nucleus), and the deviation from 
this scaling law has been precisely calculated by Karshen- 
boim [5], and more recently by Pachucki [5]. Then it is 
possible to eliminate the Lamb shift by forming a suitable 
linear combination of two frequency measurements. For 
instance, in the linear combination 7i^(2Si/2 — 8D5/2) — 
:^(lSi/2 - 2S1/2) of the 2S1/2 - 8D5/2 and IS1/2 - 2Si/2 
frequencies, the quantity Lis — 8L2S appears and the ef- 
fect of the proton charge distribution is eliminated (Lis 
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and L2S are the Lamb shifts of the IS and 2S levels). 
For example, from the measurements of the IS1/2 — 2Si/2, 
2S1/2 — 8D5/2 and 2S1/2 — I2D5/2 frequencies in hydrogen 
and deuterium, one obtains a value of Roo with a relative 
uncertainty of 9 x 10~^^ |3]. Moreover, this method gives 
the values of the Lamb shifts, and, taking into account the 
theoretical calculations of the Lamb shift, it is possible to 
deduce a value of the rms charge radius of the proton 
(rp = 0.8746 (94) fm) which is more precise than the one 
deduced from the electron-proton scattering experiments. 
In this method, the accuracy is presently limited by the 
uncertainties of the 2Si/2 — 8D5/2 and 2Si/2 — I2D5/2 fre- 
quency measurements. In our experiment ^, these uncer- 
tainties were mainly constrained by the natural width of 
8D and 12D levels and by the inhomogeneous light shift 
experienced by the atoms passing through the gaussian 
profile of the laser beams. 

To circumvent the limitation due to the uncertainties 
of the 2S— nD measurements, our group studies the IS — 3S 
two-photon transition, with the aim to deduce i?oo and rp 
from the comparison between the IS — 2S and IS — 3S 
frequencies. Indeed, as, in an atomic beam, the number 
of hydrogen atoms in the IS level is about eight orders of 
magnitude larger than the number of metastable atoms, 
the IS — 3S transition can be observed with very low light 
intensity so with a negligible light shift. In 1996, we have 
observed the IS — 3S transition and deduced, from the 
comparison with the 2S — 6S/D frequency intervals, a value 
of the IS Lamb shift J^. Since then we have undertaken 
optical frequency measurements of the IS — 3S transition. 
Because there is no simple way to determine the veloc- 
ity distribution of the IS atomic beam, a difficulty is the 
determination of the second order Doppler effect, which 
induces a red shift of —v^jlc?. To measure it, the atomic 
beam is placed in a transverse magnetic field which in- 
duces a motional electric field and so a quadratic Stark 
shift varying also as v^ . The velocity of the atoms and the 
second order Doppler effect are deduced from the vari- 
ation of this Stark shift with the magnetic field |11I12| . 
With these techniques, we have made a preliminary mea- 
surement of the IS — 3S optical frequency. 

The aim of this paper is to relate in detail these last 
measurements. Section 2 describes the experimental method. 
Section 3 and 4 are devoted to the method used to deter- 
mine the velocity distribution and to the theoretical anal- 
ysis of the line shape. Finally the results are presented and 
analyzed in section 5. 
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Fig. 1. Experimental setup for the excitation of the 1S-3S 
transition and the optical frequency measurement (Laser 
TiSa: titanium sapphire laser, FP: Fabry-Perot cavity, 
DL/Rb: laser diode stabilized to a two-photon transition 
of Rubidium). 



2 Spectroscopy of the 1S-3S transition 

2.1 Experimental setup 

The experimental setup (see Figure [T]) has been described 
elsewhere [4I10I12J . The excitation wavelength of the IS- 
3S two-photon transition is in the UV range, at 205 nm. 
This radiation is produced by quadrupling in frequency a 
CW titanium sapphire laser at 820 nm, with two frequency 
doubling stages by a lithium triborate crystal (LBO) and 
a beta-barium borate crystal (BBO). These crystals are 



placed in two successive enhancement ring cavities (LBO 
cavity and BBO cavity respectively). The first frequency 
doubling is efficient and delivers about 800 mW at 410 nm 
from an incident power of 2 W at 820 nm [l^ . This first 
cavity is locked to the laser wavelength without any length 
modulation thanks to the polarization method [13]. The 
second doubling step is by far more challenging |15) . For 
the second harmonic generation at 410 nm, the only pos- 
sible choice is a BBO crystal used at the limit of its phase- 
matching angle. This results in a low conversion efficiency, 
with a UV power below 1 mW. The second doubling step is 
operated in an O2 environment to slow the chemical reac- 
tions at the surface of the crystal. Furthermore, a photo- 
refractive effect takes place inside the crystal, which re- 
sults in a reflected blue beam at 410 nm appearing from 
one end facet of the BBO crystal. A counter propagat- 
ing wave at 410 nm develops in the ring BBO cavity. To 
reduce this effect, we have worked in a quasi-continuous 
regime where the UV intensity consists of 6 /us pulses at a 
frequency of 30 kHz. This is done by overmodulating the 
length of the BBO cavity at the frequency vq of 15 kHz, 
so as to be resonant only twice per modulation period for 
6 /zs. First-order Doppler effect due to the motion of the 
cavity mirror results in frequency shifts, which induce a 
splitting of the observed 1S-3S line. This effect is described 
in the references [12 16 . 

The frequency stabilization of the titanium sapphire 
laser is described in reference [4j. The setup is shown in 
Figure [5J The short term and long term stabilities are as- 
sured by two Fabry-Perot cavities, labelled FPA (auxiliary 
Fabry-Perot) and FP respectively. The principle of this 
stabilization arrangement is to lock the titanium-sapphire 
laser to the FPA cavity, the FPA cavity to the FP cav- 
ity and, flnally, the FP cavity to a diode laser stabilized 
on a two-photon transition of rubidium. A secondary laser 
beam from the titanium-sapphire laser is split after a dou- 
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Fig. 2. Experimental setup for the frequency stabilization 
of the titanium sapphire laser. See the explanations in the 
text (TiSa: titanium sapphire laser, AOM: acousto-optic 
modulator, EOM: electro-optic modulator, FP and FPA: 
Fabry-Perot cavities, DL/Rb: diode laser stabilized on a 
two-photon transition of Rb). 



ble pass through an acousto-optic modulator (model 3200 
from Crystal Technology at 200 MHz) and sent on the 
FPA and FP cavities. The FPA cavity (free spectral range 
600 MHz and finesse of about 400) is placed in a robust 
brass vacuum box (wall thickness of 2 cm) and carefully 
isolated from external vibrations [T3]. To reduce the fre- 
quency jitter, the laser is locked, in a first step, to the 
FPA cavity with a FM sideband method [17]. Thanks 
to this servo-loop, the frequency jitter is reduced from 
500 kHz (free running laser) to about 2 kHz |13| . 

The long term stability is guaranteed by the FP cav- 
ity. This cavity is very stable. It consists of a 50 cm long 
zerodur spacer and two silver coated mirrors, one fiat and 
one spherical (60 cm curvature radius). Its finesse is about 
120 at 800 nm. A piezoelectric transducer moves the fiat 
mirror thanks to a mechanical construction (made in fused 
silica) which avoid the rotation of the mirror (the princi- 
ple is to deform a parallelogram) [18 . This cavity is also 
placed in a vacuum box with the same design than for the 
FPA cavity. To obtain long term stability, the FP cavity 
is stabilized on a standard laser, namely a laser diode at 
778 nm stabilized to the 5Si /2-5D5/2 two-photon transi- 
tion of rubidium (DL/Rb laser). This standard has been 
described previously |19I20I21] . With this setup, as the ze- 
rodur spacer is very stable, the servo-loop of the FP cavity 
on the DL/Rb laser always uses the same fringe of the FP 
cavity and the length of the FP cavity is exactly known. 
Consequently, to excite the IS— 3S transition, the titanium 
sapphire laser is always locked on a fringe of the FP cav- 
ity which is also exactly known (fringe number 1219477 of 
the FP cavity) . The advantage of this method is that the 
1S-3S signal always appears around the same frequency 
of the acousto-optic modulator. The frequency character- 
istics of this system will be presented in section 12.31 Fi- 
nally, to scan the laser frequency, we sweep the frequency 
of the radiofrequency wave which drives the acousto-optic 
modulator. 



The 1S-3S transition is excited in a thermal IS atomic 
beam colinear with the UV laser beams. The hydrogen 
atoms are produced through a radio-frequency discharge 
from molecular hydrogen. In order to increase the UV in- 
tensity, the atomic beam is placed inside a linear build-up 
cavity (UV cavity) formed by two spherical mirrors (ra- 
dius of curvature 25 cm, 49 cm apart). Inside the cavity, 
the UV beam is focused within a waist of 48 /.tm. Because 
of the average characteristics of the UV mirrors (refiection 
and transmission of the input mirror 89 % and 8.5 %, re- 
fiection of the end mirror 96 %), the cavity finesse is about 
40. The mirrors are mounted on piezoelectric translators 
and the cavity length is locked to the 205 nm wavelength. 
This locking is detailed in reference [12]. For this servo- 
loop, the amplitude of the frequency modulation due to 
the modulation of the BBO cavity is too small with re- 
spect to the width of the UV cavity to obtain a good error 
signal. Therefore the length of the UV cavity needs to be 
modulated at the same frequency as that of the BBO cav- 
ity (15 kHz), but in quadrature. This way, the UV pulses 
issued from the BBO cavity are sent to the UV cavity 
when the length of the UV cavity is at an extremum. If 
the UV cavity length is minimum for one pulse, it will be 
maximum for the following one and so on: the UV pulses 
test successively the two sides of the Airy's peak of the UV 
cavity. For this servo loop, the modulation amplitude is a 
small fraction of the width of an Airy's peak of the UV 
cavity (typically 20 %). Finally, the UV intensity inside 
the cavity is monitored by a photodiode placed after the 
end mirror and a phase sensitive detection at 15 kHz com- 
pares the transmission of the UV cavity for two successive 
pulses to obtain the servo signal. 

The two-photon transition is detected by monitoring 
the Balmer-a fiuorescence due to the radiative decay 3S- 
2P. This fiuorescence is collected with a spherical mirror 
and a //0.5 aspheric lens system and selected with an in- 
terference filter at 656 nm. We have used two different 
detectors, a photomultiplier (Hamamatsu R943-02) or a 
CCD camera (Princeton Instruments Spec lOOB). The fol- 
lowing section presents the data acquisition in these two 
cases. 



2.2 Experimental spectra 

For the photomultiplier detection, the data acquisition 
takes benefit from the short time response of the pho- 
tomultiplier. Each scan is divided in 31 frequency points. 
For each point, the photomultiplier signal is counted dur- 
ing 1 s and we make 10 scans of the line to achieve a 
7 minutes run. Furthermore, acquisition electronics are 
designed to select only the time window during which the 
205 nm UV light is resonant inside the atomic beam cavity 
(as a consequence of the 15 kHz over modulation of the 
BBO cavity). The signal of the photomultiplier is sent, af- 
ter an amplifier discriminator and an electronic switch, to 
a multiplexer. To reduce the noise due to the dark current 
of the photomultiplier, the electronic switch transmits the 
signal only when the instantaneous UV intensity is above 
a reference level. Then, a multiplexer dispatches the signal 
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Fig. 3. Record of the 1S-3S two-photon transition in hy- 
drogen with the photomultipher (color online). This spec- 
trum is the mean of 8 runs. For each point the acquisi- 
tion time is 80 s. The two curves correspond to the sig- 
nal recorded when the BBO cavity length increases (blue 
curve, blue shifted) or decreases (red curve, red shifted). 
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Fig. 5. Schematic diagram of the optical setup for the 
optical frequency measurements (Ti:Sa: titanium-sapphire 
laser, BS: beam splitter, DM: dichroic mirror, PD: photo- 
diode). The repetition rate is detected with PDl, the off- 
set /o with PD2 and the beat note between the CW Ti:Sa 
laser and the frequency comb with PD3 (color online). 
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Fig. 4. Record of the 1S-3S two-photon transition in hy- 
drogen with a CCD camera. The total acquisition time is 
about 20 minutes. In this case the signals obtained when 
the BBO cavity length increases or decreases are not sep- 
arated. 



to two counters in phase with the modulation at 15 kHz: 
when the BBO cavity length increases, the signal is sent 
to counter 1, and, when it decreases, to counter 2. This 
way, we can observe the splitting of the IS — 3S line due 
to the Doppler shift induced by the motion of the mir- 
ror (see section [53J ■ As the signal-to- noise ratio is small, 
we take the mean of several runs to obtain an observable 
signal. Figure [3] shows the atomic signal. The two curves, 
corresponding to the two counters are separated by about 
1.1 MHz in terms of atomic frequency. Consequently we 
will have to take into account this effect in the line shape 
analysis. 



In comparison with the photomultipher, the CCD cam- 
era has two advantages: the quantum efficiency is higher 
(about 90 % for a back-illuminated camera when it is 15 
% for the photomultipher) and it is possible to obtain an 
image of the fluorescence of the atoms in the laser beam. 
The drawback is the longer time response of the CCD, 
of the order of one millisecond, which does not allow us 
to perform the electronic time selection described above. 
Moreover a readout noise, independent of exposure time, is 
also superimposed on each pixel (or group of pixels) every 
time the chip is being read. This leads to the use of long 
exposure times, of the order of 1 minute, to gather as many 
signal photons as possible before reading the chip. Never- 
theless, this exposure time should remain small compared 
to the characteristic drifting time of the titanium sapphire 
laser and of the frequency doubling stages. Finally the ex- 
posure time of each frequency point of the scan is 37 s 
and the acquisition of the 31 frequency points lasts about 
twenty minutes. Then each image is carefully analyzed to 
reduce the parasitic signal due to the UV and to take into 
account the variation of the UV intensity during the scan 
(see the references [22123) for the detail of this analysis). 
An example of IS — 3S spectrum obtained with the CCD 
camera is shown in figure 21 For this record the UV inten- 
sity was about 70 % the one of the record shown in figure 
[21 As the two-photon excitation is a quadratic process, this 
corresponds to a reduction by a factor of about 2 of the ex- 
citation probability. Nevertheless the CCD camera signal 
is about five times larger than the photomultipher signal. 
This corresponds to the gain in the quantum efficiency 
and to a better collection of the atomic fluorescence. For 
this reason, the main results presented in this article have 
been obtained by using the CCD camera. 
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2.3 Optical frequency measurement 

The optical frequency measurements are made with an 
optical frequency synthesizer, following the design intro- 
duced by Hall and Hansch J24I25) . Figure [S] shows the 
experimental arrangement. We use a second titanium sap- 
phire laser, a six chirped mirrors mode-locked femtosecond 
laser with a repetition rate frep of about 900 MHz (laser 
GigaJet 20 from Menlo Systems GmbH) with a 5 W pump 
laser (laser Verdi V5 from Coherent). The output of the 
titanium sapphire laser forms a frequency comb. The fre- 
quency /at of each line of the comb is controlled by the 
frequency rate frep and the global shift /o of the frequency 
comb with respect to the zero frequency: fpf = Nfrep + fo- 
The offset frequency /o is determined thanks to the self- 
referencing technique. The spectrum of the femtosecond 
titanium sapphire laser is broadened inside a microstruc- 
ture photonic crystal fiber (from CrystalFiber) as to span 
over more than an octave. Then the infrared part of the 
spectrum is frequency doubled in a KNbOa nonlinear crys- 
tal to obtain green radiation which is recombined, after a 
time delay, with the green part of the spectrum generated 
by the photonic crystal fiber. The result is a beat note de- 
tected by the photodiode PD2 (see figure [SJ at the offset 
frequency /o or at the frequency frep — fo- As the out- 
put spectrum of the comb before the fiber (ranging from 
780 to 820 nm) includes the wavelength of the two-photon 
excitation, the CW titanium-sapphire laser is mixed with 
the frequency comb before the photonic crystal fiber. The 
advantage is that the resulting beat note (frequency /i) is 
very stable with a good signal-to-noise ratio (40 dB in a 
300 kHz bandwidth). This frequency /i corresponds to the 
frequency difference between the titanium sapphire laser 
frequency and the closest frequency of the comb, either 
lower (case (a)) or higher (case (b)). The spacing between 
the lines of the comb is fixed by phase-locking the repeti- 
tion rate to a reference signal related to a Cs clock. Thanks 
to an optical link between our laboratory and LNE {Labo- 
ratoire National d'Essais)-SYBTE (Systeme de Reference 
Temps Espace) at the Observatoire de Paris I26j . we re- 
ceive a radio-frequency signal at 100 MHz locked to a pri- 
mary frequency standard. Then a radio frequency chain 
generates a signal at 11 GHz which is mixed with the 
twelfth harmonic of the repetition rate ^7\ . The phase er- 
ror signal is then amplified and fed back to a piezoelectric 
translator controlling the length of the femtosecond ring 
cavity. Finally the frequency of the titanium-sapphire laser 
is given by: 
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Fig. 6. Measurements of the frequency /(1219477) of the 
peak numbered 1219477 of the FP cavity, in 2005 (upper 
part) and in 2009 (lower part). The horizontal axis corre- 
sponds to the successive number of the runs (about 200 
runs during 15 days in 2005 and 300 runs during 30 days 
in 2009). Each point corresponds to a 20 minute data ac- 
quisition run in 2005, and 7 minute in 2009. During this 
period, the drift of the cavity is of a few kHz. 



standard deviation of the mean. The absolute frequency of 
the DL/Rb standard has also been measured several times. 
The results obtained after extrapolation to zero laser in- 
tensity (in order to eliminate the light shift of the two- 
photon rubidium transition) are 385 285 142 376.3 (1.0) 
kHz in May 2004, 385 285 142 376.4 (1.0) kHz in June 2005 
and 385 285 142 378.8 (1.0) kHz in February 2009. These 
values are close to the first measurement of the DL/Rb 
standa rd en 1996 which was 385 285 142 376.7 (1.0) kHz 
[20128] , 

During the recordings of the 1S-3S spectra, the fre- 
quency of the titanium-sapphire laser was continuously 
measured. Figure [6] shows the optical frequency of the 
peak of the FP cavity numbered 1219477 which is used to 
lock the titanium-sapphire laser, measured in 2005 and in 
2009. This figure illustrates the repeatability of the setup 
used to stabilized the titanium-sapphire laser (see figure 
[5]). The frequency jumps are due to some degradations of 
the servo-loops of the lasers on the FP cavity. Usually it 
is due to a bad alignment of the lasers with respect to the 
cavity or to a diminution of the electronic gain. The shift 
of about 5 kHz between the 2005 and 2009 measurements 
can be explained by the frequency of the DL/Rb laser or 
by an ageing of the silver coated mirrors of the FP cavity. 



with a -|- sign in case (a) and — sign in case (b). Then, 
by mixing both beat notes electronically and filtering, one 
can record the suitable combination of /o and /i with a 
tracking oscillator and a frequency counter. 

Thanks to the optical link with the Observatoire de 
Paris, the DL/Rb standard laser has been simultaneously 
measured during 300 s in our laboratory and in the LNE- 
SYRTE laboratory. The result of this test is a frequency 
difference of 2 (26) Hz, the uncertainty corresponds to one 



3 Determination of the velocity distribution 

In a thermal hydrogen beam at room temperature, the 
typical atomic velocity is 3 km/s. For the 1S-3S transi- 
tion, this velocity induces a second-order Doppler shift 
-i^(lS - 3S) X v^/2c^ which is about 146 kHz (i/(lS - 3S) 
is the atomic frequency). To measure this effect, we use 
a method proposed in reference |llj which is convenient 
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Fig. 7. Experimental arrangement for the determination 
of the second-order Doppler shift. The atomic beam is 
horizontal. The UV light is focused in front of the detec- 
tion system using a photomultiplier (PM). The two coils 
produce a vertical magnetic field B perpendicular to the 
direction of the atomic beam. 



when the IS-nS (or 2S-nS) two-photon hydrogen transi- 
tions are produced in an atomic beam. The principle is to 
apply a transverse magnetic field B with respect to the 
direction of the atomic beam (see figure [T]). This mag- 
netic field produces a motional electric field E = v x B 
which induces a quadratic Stark shift proportional, as the 
second-order Doppler effect, to v^ . From the observation 
of this effect, we can deduce the second order Doppler 
shift. 

The first effect of the magnetic field is a Zeeman split- 
ting (see figure E]) . In low magnetic field, because of the 
selection rules of a two-photon transition (Z\F = and 
Amp = |in]), the line between the hyperfine levels 
lSi/2(f = 1) and 3Si/2(^ = 1) is split in three lines. For 
a higher magnetic field, the line between the lSi/2(-F = 
l,mp = 0) and 3Si/2(i^ = 0,mF — 0) Zeeman sub-levels 
becomes allowed, because the hyperfine structure of the 
3S level is in the Paschen-Back regime when it is not the 
case for the IS level. Practically, the shift of the mp = ±1 
components is very small, because, in first approximation, 
the Lande factors are the same for the IS and 3S lev- 
els. Following a relativistic calculation, the Lande factors 
gj of the IS and 3S levels are respectively 2.002284 and 
2.002315 [30 . This difference induces, for a field of 20 mT, 
a residual Zeeman shift of about ±4.3 kHz. For this field 
the diamagnetic shift is about 809 Hz. On the other hand, 
the F = 1, mp = component is shifted, because the Zee- 
man effect is not the same for the lSi/2(i^ = l,mp — 0) 
and 3Si/2(F = l,mp = ) sub-levels. For 20 mT, this 
shift is about 200 MHz. 

The second effect of this magnetic field is a motional 
electric field E = v x B. For a velocity of 3 km/s and a 
magnetic field of 20 mT, this electric field is 0.6 V/cm. 
It induces a coupling of the 3Si/2 level with the 3Pi/2 




Magnetic field (mT) 

Fig. 8. Zeeman diagram of the 1S-3S transition. The 
levels are labelled by the quantum numbers [F^mp). 
In the intermediate regime the lSi/2-3Si/2(F — 1) line 
is split in four components following the selection rule 
Amp — 0. The motional Stark effect is important, for 
the lSi/2-3Si/2(i^ = 1,™^? = —1) component, around 
the crossing between the 3Si/2 {F = l,mp ~ —1) and 
3Pi/2(-F = l,m,p — 0) levels. For high magnetic field, 
the magnetic quantum numbers of these two levels are re- 
spectively {mj = —1/2, m/ = —1/2) and {mj — 1/2, 
TO/ = -1/2). 



and 3P3/2 levels. For the 3Si/2 level, in the limit of high 
magnetic field, the electronic and nuclear spins are de- 
coupled and the good quantum numbers are (J, /, m.j, 
m,i). For example, the 3Si/2(-F = l,mp — —1) level 
corresponds to the quantum numbers mj — —1/2 and 
TO,/ = —1/2. It is mainly coupled to the 3Pi/2(to.j — 
1/2, mi = -1/2), 3P3/2(to,/ = 1/2, to/ = -1/2) and 
3P3/2('TIj = —3/2, m/ — —1/2) sub-levels following the 
selection rules Amj = ±1 and Ami = 0. This kind of mo- 
tional Stark mixing is well known and, for instance, was 
used by Lamb and Retherford to polarize a metastable 
atomic beam ^HT. We consider the quadratic Stark ef- 
fect due to this motional electric field. This shift, propor- 
tional to v^, is very small, but observable in the case of 
the mp = ±1 line components, because the Zeeman shift 
of these components is negligible. In particular, this ef- 
fect is important around the level crossing between the 
3Si/2(F = 1,tof = -1) and 3Pi/2(F = 1,tof = 0) 
Zeeman sub-levels, which appears for a magnetic field of 
about 18 mT (see figure |S]). 
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The complete calculation of the line shape will be de- 
scribed in the section |4l We present here a simple picture. 
The initial state g (the levels lSi/2(^ — l,mp — ±1)) 
is coupled by two-photon excitation to an excited state e 
(here the levels 3Si/2(-F' = 1,tof = ±1))- This state e is 
mixed by the Stark hamiltonian Vef with several / states 
(the Zeeman sub-levels of the 3Pi/2 and 3P3/2 states). In a 
simple approach, the perturbed energy E'^ and width F^ of 
the considered e level are deduced from the unperturbed 
energy and width E^ and Fe, using the perturbation the- 
ory in lowest order: 



E'^-^^^^. = E,.^^^^,. + Y: 



Wr- 



ef\ 



[E,-Ef)-i'^{F,-Ff) 



Then the Stark shift 5s of the level e is given by: 



(2) 



5s = Re(S^ - E, 



E 



{E,~Ef)\V,j\ 



{E,~EfY + {%Y{F,-FfY 



In these equations, the energies E^ and Ef take into ac- 
count the shifts of the e and / levels due to the Zeeman 
effect: then £^e — Ef \s a. function of the magnetic field 
and the Stark shift appears as a sum of dispersion curves 
corresponding to the different level crossings between the 
3Si/2 and 3Pj levels. There is no divergence at the level 
crossing when E^ — Ef because of the difference between 
the natural widths of the 3Si/2 and 3Pi/2 levels (respec- 
tively 1 MHz and 30.6 MHz). As the effect of the electric 
field is negligible for the IS1/2 level, the total shift 5 of 
the two-photon line is: 



5 = 5z + 5s- i^(lS - 3S) X 



2c2 



(4) 



where 5z is the shift due to the difference of the Zeeman 
effect for the g and e levels and the third term the second 
Doppler effect for an atom of velocity v. 

Figure IHl shows, for an atom at 3 km/s, the shifts of 
the lSi/2— 3Si/2(i^ — l,mF = ±1) transitions around the 
level crossing at 18 mT. The large dashed line represents 
the line position without any shifts for the velocity v = 
and the small dashed line the line position shifted by the 
second-order Doppler effect for u = 3 km/s and B = 0. 
This shift is about —146 kHz. With respect to this refer- 
ence, the shift of the lSi/2-3Si/2(-F = l,mi? = -1-1) line 
(curve (a)) is small (about 10 kHz), because, in this range 
of magnetic field, the 3Si/2(-F' = l,mF = +1) level is 
far from the 3Pi/2 and 3P3/2 levels. On the other hand, 
the shift of the lSi/2-3Si/2(i^ = l,mi? — —1) line (curve 
(b)) is important around the level crossing. It has a dis- 
persion shape and is larger than the second-order Doppler 
effect. Moreover, there are two values of the magnetic field 
(points A and B) where the second-order Doppler effect 
is exactly compensated by this motional Stark effect. For 
these values, the position of the line is independent of the 
atomic velocity, because the second-order Doppler effect 
and the motional Stark effect both vary as w^. Neverthe- 
less, this effect cannot be used directly, because the IS1/2- 
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Magnetic field (mT) 

Fig. 9. Shift of the 1S-3S line due to the motional elec- 
tric field (calculation for an atom at a velocity of 3 km/s) . 
Zero: line position for an atom at rest; small dashed line: 
shift of the line due to the second-order Doppler effect. 
Curves (a) and (b): line position of the lSi/2-3Si/2(F = 
l,mF = +1) and lSi/2-3Si/2(^ — I^ttif = —1) com- 
ponents. In A and B, the second-order Doppler effect 
of the lSi/2-3Si/2(^ = 1,tof = —1) component is ex- 
actly compensated by the motional Stark effect. Curve 
(c): global motional Stark effect for the two components 
lSi/2-3Si/2(i^ = l,mf = ±1). 



3Si/2(-F — \,niF — ±1) lines are not separated (the mo- 
tional Stark effect is small with respect to the 3Si/2 natu- 
ral width (about 1 MHz). Curve (c) shows the barycenter 
of these two lines. Though reduced by a factor two, the 
dispersion amplitude is about 230 kHz, and the compen- 
sation of the second-order Doppler effect is still 73% for a 
magnetic field of 17 mT. We have observed this effect to 
deduce the second-order Doppler effect. 



4 Line shape 

4.1 Theoretical background 

The aim of this calculation is to take into account simul- 
taneously the natural width and the Zeeman and Stark 
effects. We follow the procedure described in reference [4]. 
The evolution of the density operator p is: 

'^ = \[{H^ + Vl + Vs),p\+Fp (5) 

at in 

where Vs is the Stark Hamiltonian and the operators Vl 
and F describe the two-photon excitation and the spon- 
taneous emission. We use the notations of section [3] the 
states g, e and / are respectively the IS1/2, 3Si/2 and 3Pj 
levels. We make the rotating wave approximation and we 
introduce the two-photon Rabi frequency fie'- 



(e| Vl Iff) 



TiQp 



eyi^{—2iujt) 
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n. 



8™g|(e|Qtp|<?)|/ 



(6) 



where w/27r is the laser frequency, / the power density of 
the Ught, flo the Bohr radius, a the fine structure constant 
and m the electron mass. For a polarization along the z- 
axis, the two-photon operator Qt^ is given in atomic units 
(?i = ac = TO = 1) by: 



^ — ^ ^J — I 



(7) 



urg 



where the sum is made on all the atomic states r and 
tjJrgl^'K is the atomic frequency difference between the g 
and r states. Because of the selection rules for the two- 
photon excitation, the initial state g is coupled to a single 
excited state e. 

If we assume that J7e ^ ^e, we can neglect in a first 
step the populations and coherences pee', P//' or p^f of 
the upper levels. In the rotating frame, we replace the 
density operator by an operator a with Ogg — pgg, a eg — 
Peg exp(2ia;i) and age = Pge exp(— 2ia;t) and we introduce 
the frequency detunings Ae = 2uj — (we — Wg) and Af — 
2uj — [ujf ~ ujg) {huji is the energy of the level i). In this 
way, we obtain from equation ([5]) a set of equations: 



d(7n 



dt 



da, 



eg 



dt 



Fe 



.a 



.) 



(8) 



iAe-^j aeg - i^dgg - T XI ^e/^/9 (9) 






(10) 



where Vfe = {f\Vs\e). Then we assume that the optical 
coherences aeg and afg follow adiabatically the population 



agg, i.e. that: 

^ = and ^ = 
dt dt 

With these hypotheses, the equations ^ [TU)) give: 



iAe-^\ aeg + X 



/ 



0^/ - ^) 



Oen = J-I-Cr 



h'\iA 



'eg 



99 



This equation gives the coherence aeg as a function of the 
population agg. It can be written more simply: 



where: 



A^XiAe- 



Fe 



1 .^e 



E' 



VefVfe 



2/ yn^ 



0^/ - ^ 



(12) 



(13) 



Combining with equation ([8]) , one obtains the evolution of 
the population agg and the probability of the two-photon 
excitation. 



To obtain the correct line shape, we have to calculate 
the 656 nm fluorescence from the n — Z levels towards the 
n — 2 levels. This fluorescence is due to the cascades 3Si/2 
-^ 2Pj and 3Pj -^ 2Si/2. Because of the Stark mixing, 
the 3Pj states are also populated. Then, to obtain the line 
shape, it is required to calculate the populations aee. Off 
and affi . They are deduced from equation [5j 



da„ 



, = —TeOee + ^ ^e {^eg ~ Oge) 



/ 



(TefVfe) 



(14) 



doje 

dt 



re + Tf 



T^eOfg 



-LOe)crfe 
Vfe (o-ee " Off) 



afe - I [UJf 

i 



da 



ff 



dt 



-^f^ff ~ ^ (VfeOef - OfeVef) 



(15) 



(16) 



In our experimental conditions, the two-photon excitation 
probability is very small, typically lO^'^/s for an atom at 
the center of a 1 mW UV beam. As this probability is very 
small with respect to the natural width of the 3S level (1 
MHz), we are in a stationary regime and daw/dt = 
(with i and i' are the levels e or /). Then it is possible 
to calculate the populations aai from the set of equations 

dH-nni). 

First we consider the simple case where the e level is 
coupled with only one / level by the Stark Hamiltonian. 
From equations (jlOp . (TT^ and (TTSl) . wc deduce: 



J-Vfe {Oee " Off) 



Vf. 



Ofe 



' ih A Af+irf/2^ 



UJf 



;£i±£L 



(17) 



The populations aee and Off are obtained from the equa- 
tions dni) and (IT6t: 



FeOe 



2 Re 



f2e 

i — ap 
2 



Re (iVefOfe) 



FfOff = -ReiiVefOfe) 



(18) 



Using the expressions of (Teg (equation (J12p ) and Ofe (equa- 
tion (J17p ). we can express aee and Off as a function of Oggi 



^ B 



1 -1 
-1 1 



Oee 



where the coefficients B, Ci and C2 are: 

2,iVef\^^re + rf) 



C2 



'99 



(19) 




C2 = ^Re 
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Finally, the populations a^e and ct// are: 



fee 



Fe- B B 
B Fj-B 



C2 



'aa 



(21) 



The calculations are similar in the case where there are 
several / levels mixed by the electric field with the e level. 
They are described in detail in the reference [5^ . 



4.2 Theoretical line shape 

For an atom in the initial state g with the velocity v in the 
magnetic field B, the detected fluorescence Fg{uj,v,B) is: 



Fg{LO,V,B)^D Te^e 



f I 



\ 



(22) 



where D describes the efficiency of the photon detection 
and 7/ is the branching ratio of the fluorescence from the 
3P levels towards the 2S levels (about 0.1183). The popu- 
lations are calculated following equation (|2T|) and, to take 
into account the second-order Doppler effect, the laser fre- 
quency u is replaced by a;(l -|-w^/2c^) in the expression of 
Z\e and Z\/. FigurefTUlshows the line shape of the two com- 
ponents rap = 1 and rtip = — 1 for a velocity of 3 km/s and 
a magnetic fleld of 17.1 mT. This calculation takes into 
account the Stark mixing of the 3Si/2(-F = l,mp = ±1) 
with all the 3Pi/2 and 3P3/2 levels. For the mp = 1 com- 



ponent, the effect of the Stark mixing is negligible and the 
line is red shifted by the second-order Doppler effect. On 
the other hand, the mp = — 1 component is blue shifted 
and the shift due to the Stark mixing is larger than the 
second-order Doppler effect. The widths of the two com- 
ponents are also different: 1.451 MHz for the mp = —1 
component and 1.006 MHz for the mp — 1 one (the nat- 
ural width of the 3S level is 1.005 MHz). 

In the precedent calculation, we have considered an 
atom which is continuously inside the UV beam and equa- 
tion (im gives the probability per unit of time that we de- 
tect a 656 nm photon. Now we have to take into account 
the velocity distribution of the atoms. The atomic veloc- 
ity distribution of a hydrogen beam produced by a radio 
frequency discharge dissociator has been studied by Jadus- 
zliwer and Chan [33 . In our experimental conditions (dis- 
sociator pressure in the range of 0.4 torr), we can assume 
that this distribution f{v, a) is close to the Maxwellian 
form: 

f{v, a) = v^ exp{~v^/2a^) (23) 



where a — yJkT /M (T temperature, M atomic mass) 
Then the line shape R{i^, a, B) is: 



R{uj. 



a,B)^a -f{v,a)y2Fg{uj,v,B)dv 

JO ^ 



(24) 



where g are the levels 3Si/2(^ = l,mp = ±1) and a a 
normalization factor which depends of the atomic density 
in the atomic beam. If L is the length of the part of the 
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Detuning from the atomic frequency (IVIHz) 

Fig. 10. Line shape of the mp — +1 and m,p ~ —1 Zee- 
man components of the line for an atom at a velocity of 
3 km/s in a magnetic field of 17.1 mT. The dashed line 
is the line position for an atom at rest. The red shift of 
the mp = +1 component is principally due to the second- 
order Doppler effect. On the contrary, because of the mo- 
tional Stark effect, the m,p = — 1 component is blue shifted 
and broadened and its intensity is reduced. 



atomic beam in front of the detection system, the two- 
photon excitation probability is proportional to the transit 
time L/v. Consequently, we have introduced the factor 
1/v into equation ([M)) . In our calculation, we have also 
taken into account a slight dependence of the two-photon 
excitation probability with the velocity v which is due to 
the form of the UV gaussian beam and the life time r of 
the excited level: if the detection is made exactly at the 
waist of the UV beam, the excitation has been made at 
a distance vt from this waist. Then, in equation ([B]), the 
power density of the light / depends on the velocity v. 
The calculation of this effect is described in detail in the 
reference [23] . 

As explained in section 12.11 the length of the BBO 
cavity is overmodulated: the UV intensity consists in a 
succession of pulses which corresponds to an increase or a 
decrease in the cavity length. A splitting and a broadening 
of the line consequently appear. This effect is described in 
reference |T6^. Following [34], the line shape of the two- 
photon transition is given by: 



+00 



F{Q)\'R{^,a,B]dn (25) 



where F{Q) is the Fourier transform of the square of the 
electric field at 205 nm. 

If ujlI'^t^ is the laser frequency at 410 nm inside the 
BBO cavity, the electric fleld is: 



Eiio(t) = Re 



TcEoex.p{iujLt) 



1 — re expiiSt) 



(26) 



where Sqc*"^' is the incident fleld on the BBO cavity, 
Tc the transmission of the input mirror, re the reflection 
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B=17.1 mT 
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Detuning from the atomic frequency (IVIHz) 

Fig. 11. Theoretical line shape of the lSi/2(-F' — l,mp — 
±l)-3Si/2(-F — l,mp = ±1) transition for a thermal 
atomic beam (a — 1.6 km/s) in a magnetic field of 17.1 
mT. The broadening and the shift due to the modulation 
of the length of the BBO cavity are taken into account 
with the parameter A — 57 kHz. The full line is the sum 
of the two contributions due to an increase or a decrease 
in the cavity length (dashed lines). 



coefficient for a round trip in the cavity and S/2t: the 
frequency shift due to the motion of the mirror {5 > 
if the cavity length decreases) . The equation ([25)1 is valid 
because the lifetime of a photon inside the BBO cavity 
(about 0.2 lis) is short in comparison with the duration 
of the UV pulse (6 fj,s). Then the electric field at 205 nm 
generated by the BBO crystal is proportional to: 



£^205 (t) - Re 



TcEoexp{iujLt) 



I — re exp{i6t) 



Then the function F{f]) is: 



F{f2) 



1 

2^ 



T^£'(*exp[i(J7-4a;L)i] 
4 [1 — re exp{—iSt)] 



dt. 



(27) 



(28) 



where we have only kept the term resonant with the two- 
photon transition. As the duration of the UV pulses (6 //s) 
is short with respect to the period of the modulation (67 
/zs) , we can linearize the denominator in equation (j28|) and 
obtain \F{fi)\ . When the length of the cavity decreases 
((5 > 0), we have: 



\F{f2f 



if 

( rhEt 
\ 24r4 54 



c" 



Q < 4:UJL, 



in-AwL) 



X exp 



Q - Aujl 



if n> 4wl(29) 
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where A — re (5/2(1 — re) characterizes the shift and 
the broadening of the line. When the cavity length in- 
creases, the function |_F(i7)| is symmetric with respect dition to Z\ and z^eg , the other parameters of the fit are the 
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Atomic frequency - 2 922 742 930 (in IVIHz) 

Fig. 12. Fit of the theoretical curve to the experimental 
data of figure |4l The velocity distribution width is fixed 
at a = 1.6 km/s. The fitted parameters are a modulation 
A of 69 kHz and a frequency of the 1S-3S transition of 
2 922 742 936.745(17) MHz. 



to f2 = Alul. This symmetry supposes that the frequency 
shifts Sincr and Sdecr are the same for the forward and 
backward scanning. In fact, a non-linear refractive index 
of the BBO crystal can produce a dissymmetry. This ef- 
fect will be estimated in section [5] From equations (1^^ . a 
straightforward calculation shows that the mean position 
of each function |-F(^)| is shifted by ±7Z\ with respect 
to 4WL IS2]- 

Finally, the line shape is obtained from equations (P5|) 
and (P^ . Figure [11] shows the profile calculated with a 
velocity distribution parameter ct = 1.6 km/s, in a mag- 
netic field of 17.1 mT and for a modulation Z\ = 57 kHz. 
This modulation splits the line shape in two contributions 
corresponding to an increase or a decrease in the cavity 
length (dashed lines on figure [TT|) with a separation of 
about 800 kHz. 



4.3 Fit of the experimental profiles 

The parameters of the theoretical curves are the frequency 
of the 1S-3S transition ly^g — (we — "^g)/^ for an atom at 
rest without magnetic field, the parameter A which char- 
acterizes the splitting due to the modulation of the length 
of the BBO cavity, and a which describes the velocity 
distribution. In a simple model the second-order Doppler 
shift is given by — 3/2(cr/c)^i^eg g]. Consequently the pa- 
rameters i>eg and a are strongly correlated and it is not 
possible to deduce these two parameters from the fit of the 
line shape. Practically, a is fixed and the fit gives A and 
i^eg] u will be deduced in the next section from the com- 
parison of the results for different values of the magnetic 
field. 

Figure [12] shows an example of the fit of the theoretical 
curve to the data obtained with the CCD camera. In ad- 
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Fig. 13. Fit of the theoretical curves to the experimental 
data of figure [31 The upper part of the figure (respectively 
the lower part) shows the fit of the sum (respectively the 
difference) of the two curves of figure[3]corresponding to an 
increase or a decrease in the length of the BBO cavity. The 
velocity distribution is fixed at cr = 1.6 kni/s. The fitted 
parameters are a modulation A of 92 kHz and a frequency 
of the 1S-3S transition of 2 922 742 936.722(18) MHz. 



amplitude of the signal and an offset corresponding to var- 
ious detection noises (parasitic light and read-out noise of 
the CCD camera) . The signal-to- noise ratio of this record 
is about 24 and the atomic frequency is determined with 
a statistical uncertainty of about 17 kHz. 

Through detection with the photomultiplier, it is pos- 
sible to fit simultaneously the two spectra corresponding 
to an increase or a decrease in the BBO cavity length. Nev- 
ertheless it is preferable to fit the sum and the difference 
of the two curves, because, in the case of the difference, 
a large part of the noise due to the UV parasitic light is 
rejected. An example is shown on figure [T51 The signal-to- 
noise ratio is about 16 for the fit of the sum of the curves 
and 13 for the difference. The statistical uncertainty of the 
atomic frequency v^g is 18 kHz. 



5 Results 

The main results have been obtained by using the CCD 
camera. The two-photon transition has been observed for 
four values of the magnetic field, 0.029 mT, 16.01 mT, 
17.12 mT and 19.15 mT. The first value corresponds to 
the residual earth magnetic field. The magnetic field is cal- 
ibrated with an accuracy of about 0.01 mT by observing 
the Zeeman effect of the lSi/2-3Si/2(-F' = 1, f^F — 0) tran- 
sition. Table [T] gives, for each value of the magnetic field, 
the number of recordings n, the mean v^g of the atomic 
frequencies obtained by fitting each recording with a = 
1.6 km/s and the values of x^/in-l)- For B = 0.029 mT, 
the uncertainty of the frequency measurement is 5 kHz, 
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Fig. 14. Frequencies I'eg as a function of the parameter 
a for four values of the magnetic field. The calculation 
has been made for 6 values of a spanning between 1 km/s 
and 2 km/s (small dots). The full lines are the interpola- 
tions between these points and the dashed lines the error 
bars. The large dot is the best value for a and Veg with 
corresponding uncertainties (color online). 



i.e. a relative uncertainty of 1.7 x 10~^^. Unfortunately it 
is not the final accuracy, because the parameter a of the 
velocity distribution has been fixed to 1.6 km/s. 

To determine a the mean frequency v^g is calculated 
for different values of a (in the range 1 to 2 km/s) to 
obtain four curves Vf.g{<J,B) (see figure [T4l) . The values 
of (T and Veg are obtained from the crossing of these four 
curves. With a least squares method, we obtain: 

cr = 1.646 (89) km/s (30) 

Veg = 2 922 742 936.7275 (120) MHz (31) 

The corresponding point is indicated in figure [T3] (large 
dot); Vf.g is the frequency of the lSi/2-3Si/2(-F — 1) two- 
photon transition for an atom at rest. The statistical un- 
certainty (12 kHz) corresponds to a relative uncertainty 
of 4.1x10-12. 

For a thermal atomic beam one has a — ^JkT/M (see 
section 1121) ■ The value of a (equation (UHl)) corresponds 
to a temperature of 328(35) K which is compatible with 
a hydrogen atomic beam produced with a radio frequency 
discharge. This a value is also in agreement with our pre- 
vious measurements |12] . In this work, we did not mea- 
sure the absolute frequency of the 1S-3S transition. Nev- 
ertheless, by comparing measurements with and without 
magnetic field, we had observed the effect of the motional 
electric field and deduced a value for a of 1.55(11) km/s. 
For this analysis, we had approximated the theoretical 
line shape by a Lorentzian shape. We have remade a com- 
plete analysis of these data following the procedure de- 
scribed in sectional The result is cr = 1.633(180) km/s. 
This value is more reliable than the result given in ref- 
erence |12| and in good agreement with the new result 
(equation ((50)) ). It is also possible to compare this result 
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Table 1. Features of the recordings for the different values of the 


magnetic field. 


Magnetic field 0.029 mT 


16.01 mT 


17.12 mT 19.15 mT 


Number of recordings 39 
Mean frequency v^g - 2 922 742 930 (MHz) 6.724(5) 
xV(n - 1) 0.69 


19 

6.728(10) 

1.03 


66 28 

6.721(6) 6.707(10) 

1.21 0.56 



Table 2. Error budget. 



Frequency measurements 


8 X 0.3 kHz 


Light shift 


0.3 kHz 


Pressure shift 


1.2 kHz 


Velocity distribution 


3.0 kHz 


Scan of the BBO cavity 


2.6 kHz 


Statistic 


12.0 kHz 


Quadratic sum 


13.0 kHz 



with the measurement of the velocity distribution of the 
metastable 2S atomic beam used to observe the 2S-nS and 
2S-nD two-photon transition in hydrogen [3] . In these ex- 
periments, the metastable atomic beam was obtained by 
electronic excitation of a IS atomic beam which had the 
same design than the atomic beam used in the present 
experiment. Because of the electron impact, there was an 
angle 9 = 20° between the IS and the 2S atomic beams 
and the metastable atoms were slowed with respect to 
the IS atomic beam by a factor of about cos{6). Using 
Doppler spectroscopy of the Balmer-a line, the a value 
of the metastable atomic beam was 1.525(10) km/s [55] . 
Taking into account the factor cos{6), this corresponds to 
a value of 1.623(11) km/s. In spite of the roughness of 
this model, this value is also in perfect agreement with 
the present result. 

The error budget is summarized in Table [5] The fre- 
quency measurements are made with respect to the 100 MHz 
reference signal from LNE-SYRTE. During the measure- 
ments of the 1S-3S transition, this frequency reference was 
obtained from a hydrogen maser numbered 40805 of LNE- 
SYRTE. Simultaneously to the 1S-3S measurement this 
standard was also compared to UTC-OP (Coordinated 
Universal Time-Observatoire de Paris) and linked to the 
SI (Systeme International). The result was: 



where /3g and /3e are the matrix elements of the light shift 
operators in atomic units for the level g and e. Their values 
are: /3is = -6.44539 and /^as = 20.9264 ^. The UV 
power inside the atomic beam cavity is estimated at the 
maximum to 2 mW. This produces a light shift of 280 Hz, 
rounded to 300 Hz in table [5] 

The pressure in the vacuum apparatus is about 6 x 
10~^ torr. The collision shift of the 1S-3S transition is not 
known. Nevertheless it would be similar to the ones of the 
2S-3P hydrogen transition because these shifts are domi- 
nated by the two upper levels of these transitions of which 
the wave functions have the same spatial size. The shift of 
the 2S-3P transition in He buffer gas has been measured 
by Weber et al. to be about -9 MHz/torr ^37]. Then, to 
estimate the uncertainty due to the pressure shift, we sup- 
pose an upper limit of 20 MHz/torr, corresponding to the 
value of 1.2 kHz in tabled 

Several other effects can modify the frequency deter- 
mination. In the theoretical analysis, we suppose that the 
two Zeeman sub levels lSi/2(F — 1, nip = ±1) are equally 
populated. Actually, for a magnetic field of 20 mT, the 
energy splitting between these two levels is about 10"'' 
times the thermal energy. That produces a small pop- 
ulation difference and a frequency shift of about 4 Hz 
which is negligible. To test the sensitivity of the result 
to the form of the velocity distribution, an analysis of the 
data has been made with a different velocity distribution 



^SYRTE = vsi [1 + (575 ± 6) X 10-15] 



(32) 



That produces a correction of -f 1680(18) Hz of the IS- 
3S frequency. On the other hand, the Allan variance of 
the frequency measurement of the titanium sapphire laser 
(about 5 X 10"'^'^ in 1000 s "23]) corresponds to a statistical 
uncertainty of 180 Hz for a frequency measurement during 
a 20 minute run, conservatively rounded to 300 Hz in table 

m 



Following the notations of the reference 
shift Avis is: 



Avis = 



2al 



We - Pa) , 



the light 



(33) 



f{v,a) 



■exp(-w^/2cr^). The result of equation (I5T1) 



is then shifted by —3 kHz. We have adopted this value 
as the upper limit for the uncertainty due to the veloc- 
ity distribution. A last possible effect is a dissymmetry in 
the modulation of the BBO cavity, if the velocity of the 
mirror is not the same when the cavity length decreases 
or increases. The two Doppler shifts 5decr and 5incr due 
to the motion of the mirror cavity would not be exactly 
opposite. The maximum relative difference between these 
shifts is estimated to 2%. The analysis of the data with 
this dissymmetry shifts the frequency of the 1S-3S transi- 
tion by 2.6 kHz. This value corresponds to the uncertainty 
labeled "Scan of the BBO cavity" in Table H 

Taking into account the statistical uncertainty (12 kHz), 
the final result is: 

V [IS - 3S(i^ = 1)] = 2 922 742 936.729 (13) MHz (34) 

The total uncertainty is 13 kHz, i.e. a relative uncertainty 
of 4.5 X 10-1^. After the measurement of the 1S-2S tran- 
sition, this result is the most precise value of an optical 
frequency in hydrogen. The hyperfine structure of the IS 
level is well known [38] and that of the 3S level is evalu- 
ated from the Fermi formula with the Breit correction to 
be 52 609.4 kHz. After correction of the hyperfine split- 
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ting, the frequency of the 1S-3S transition is: 

V (IS - 3S) = 2 922 743 278.678 (13) MHz (35) 

Unfortunately, the accuracy of this result is not suffi- 
cient to improve the determination of the Rydberg con- 
stant. Several methods can be used to extract R^o- The 
first one is to use the measurement of the charge distri- 
bution of the proton from electron scattering to calculate 
the Lamb shifts of the IS and 3S levels. One then obtains 
for the Rydberg constant: 



i?oo = 10 973 731.568 75 (18) m" 



(36) 



For the calculation of the Lamb shift, we have taken into 
account all the terms given in the CODATA report [5]. 
The fundamental constants used are also the values of CO- 
DATA, except for the fine structure constant. For that we 
use the recent result of Gabrielse [SS] • With a relative un- 
certainty of 1.7 X 10~^^, this result is in acceptable agree- 
ment with the value given by the CODATA adjustment of 
the fundamental constants: i?oo = 10 973 731.568 527 (73) m 

As explained in section (TJ another method to obtain 
the Rydberg constant is to compare the 1S-2S and 1S-3S 
frequencies by using the scaling law of the Lamb shifts 
[515] . The result is: 



i?„o = 10 973 731.568 



(37) 



The relative uncertainty of this result is 6.2 x 10^^^ and it 
is about ten times less precise than the CODATA value. 
The reason is that the relative accuracy of the 2S-3S tran- 
sition which appears in this calculation is only 2.8 x 10^^^. 
With this method, one also obtains a value for the radius 
of the charge distribution of the proton: 



0.911 (65) fm 



(38) 



This value is in agreement with the result deduced from 
the electron scattering (rp — 0.895 (18) fm) and from 
the adjustment of the data in hydrogen and deuterium 
(rp = 0.8760 (78) fm). 



6 Conclusion 

To conclude, the optical frequency of the 1S-3S transition 
has been measured for the first time by using a femtosec- 
ond frequency comb with a relative accuracy of4.5xl0~^^. 
It is the best measurement of an optical frequency in hy- 
drogen after the one of the 1S-2S transition. The second- 
order Doppler effect has been determined from the ob- 
servation of the motional Stark effect due to a transverse 
magnetic field. A careful theoretical analysis has been pre- 
sented to describe the main features of the line shape. In 
spite of this high accuracy, this result does not improve 
significantly the determination of the Rydberg constant. 
For that, an accuracy of few 10~^^ would be useful. 

Presently the experiment is mainly limited by the low 
intensity of the UV source at 205 nm. To circumvent this 
difficulty we plan to modify the UV source by replacing the 



two frequency doubling stages at 820 nm and 410 nm by 
a frequency sum of a UV source at 266 nm and a titanium 
sapphire laser at 896 nm. With this scheme, similar to 
the one used by Bergquist to obtain 194 nm [30], the UV 
power would be increased by a factor of ten and the signal- 
to-noise significantly improved. Moreover it would be also 
possible to observe the 1S-4S transition in hydrogen which 
lies at 194.5 nm. 

The authors thank G. Hagel for his essential contribu- 
tion to the first step of this experiment and P. Clade for 
fruitful discussions on the manuscript. They thank also O. 
Acef for the frequency measurements of DL/Rb standard 
and they are indebted to the Service du Temps of LNE- 
SYRTE laboratory for the time reference. This work was 
partially supported by the Bureau National de Metrologie 
(now Laboratoire National d'Essais). 
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